Nonlinear Reduced Order Random Response Analysis of 
Structures with Shallow Curvature 


Adam Przekop 1 

National Institute of Aerospace, Hampton, VA 23666 
and 

Stephen A. Rizzi 2 

NASA Langley Research Center, Hampton, VA 23681 

The goal of this investigation is to further develop nonlinear modal numerical simulation 
methods for application to geometrically nonlinear response of structures with shallow curvature 
under random loadings. For reduced order analysis, the modal basis selection must be capable of 
reflecting the coupling in both the linear and nonlinear stiffness. For the symmetric shallow arch 
under consideration, four categories of modal basis functions are defined. Those having symmetric 
transverse displacements (ST modes) can be designated as transverse dominated (ST-T) modes and 
in-plane dominated (ST-I) modes. Those having anti-symmetric transverse displacements (AT 
modes) can similarly be designated as transverse dominated (AT-T) modes and in-plane dominated 
(AT-I) modes. The response of an aluminum arch under a uniformly distributed transverse random 
loading is investigated. Results from nonlinear modal simulations made using various modal bases 
are compared with those obtained from a numerical simulation in physical degrees-of-freedom. 
While inclusion of ST-T modes is important for all response regimes, it is found that the ST-I modes 
become increasingly important in the nonlinear response regime, and that AT-T and AT-I modes 
are critical in the autoparametric regime. 
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Nomenclature 

mass matrix (physical and modal coordinates) 
proportional damping matrix (physical and modal coordinates) 
linear stiffness matrix (physical and modal coordinates) 
displacement response vector (physical and modal coordinates) 

force excitation vector (physical and modal coordinates) 
nonlinear restoring force (physical and modal coordinates) 
modal basis function matrix 

horizontal and vertical displacement, and rotation (global physical coordinates) 
time 

linear, quadratic nonlinear, and cubic nonlinear modal stiffness coefficients 

identity matrix 
viscous damping factor 
undamped natural frequencies 


I. Introduction 

j ^ IRECT numerical simulation simulation of nonlinear random response in physical degrees-of-freedom (DoFs) 
is computationally intensive for even the simplest structures. Its use for design of high-cycle-fatigue tolerant 


aerospace vehicle structures is considered impractical. Accordingly, much effort has been spent in recent years to 
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develop accurate reduced order analyses, such as finite-element-based nonlinear modal numerical simulation, which 
could be suitable for use in design environments. 

Reduced order finite element methods may be viewed as being in one of two categories; those in which the 
nonlinear modal stiffness is directly evaluated from the nonlinear finite element stiffness matrix (so-called direct 
methods), and those in which the nonlinear modal stiffness is indirectly evaluated. Direct methods are typically 
implemented in special purpose finite element codes in which the nonlinear stiffness is known . 1 ' 4 The only known 
implementation of a direct method in a commercial code is due to Bathe and Gracewski . 5 Indirect stiffness 
evaluation methods are typically implemented for use with commercial finite element codes in which the nonlinear 
stiffness is unavailable . 6,7 For both direct and indirect stiffness evaluation approaches, the accuracy of the solution 
depends on the selection of the modal basis, through which the nonlinear modal stiffness may be determined. The 
majority of research performed in the field of reduced order numerical simulation of thin-walled aerospace 
structures has been primarily focused on isotropic planar configurations . 1,6 " 12 Through comparison with numerical 
simulation in physical DoFs, the authors recently demonstrated the ability of a reduced order method to accurately 
predict geometrically nonlinear random displacement and stress response for a planar isotropic structures, provided 
that a suitable modal basis is utilized . 8 However, most practical aerospace applications involve more complex 
geometries. Thus, the focus of this study is to extend the authors’ previous work using an indirect stiffness 
evaluation approach to determine the effect of modal basis selection on the response of structures with shallow 
curvature. It is known that such structures may exhibit not only a hardening but also a softening type of 
nonlinearity 1 ’ 14 which may be non-conservative relative to a linear solution. 

For planar configurations, the linear bending and membrane responses are uncoupled. Only the large deflection 
nonlinear response couples the bending and membrane behavior. Therefore, basis functions derived from a linear 
eigenanalysis are readily separable into eigenvectors associated with only the bending response and those associated 
with only the membrane response. Curved structures, on the other hand, exhibit an inherent linear coupling between 
transverse and in-plane displacements. This coupling is independent of the response magnitude. Similar to the 
planar configuration, there also exists a coupling of transverse and in-plane displacements resulting from large 
deflection nonlinearity. Therefore, for reduced order analysis, the modal basis selection must be capable of 
reflecting the coupling in both the linear and nonlinear stiffness. For the structures under consideration, two 
categories of modes are defined. Those having symmetric transverse displacements, with anti-symmetric rotations 
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and in-plane displacements, are subsequently referred to as ST modes. Those having anti-symmetric transverse 
displacements, with symmetric rotations and in-plane displacements, are referred to as AT modes. The ST and AT 
modes may be further classified according to their dominant behaviors. Mass-normalized eigenvectors dominated 
by the transverse component will be referred to as a type T mode and those dominated by the in-plane component 
will be referred to as a type I mode. Consequently, this study is limited to shallow curvatures, allowing the 
identification of basis functions that are either type T or type I. For the structure considered, type T modes typically 
have lower natural frequencies than type I modes. Flence, four classes of modes are considered for inclusion in the 
modal basis, i.e. ST-T, ST-I, AT-T, and AT-I. 

The presence of non-linear coupling also implies a possibility of autoparametric resonance, 15 which is also 
referred to as autoparametric interaction or internal resonance. 16 An external loading can directly excite only a 
motion which is consistent with its characteristics, i.e., symmetric loading acting on a symmetric structure directly 
excites a symmetric transverse and anti-symmetric in-plane response. Flowever, if sufficiently high nonlinear modal 
interactions exist, a response behavior not conforming to an external loading may also be initiated. 17 In other words, 
an energy flow can occur from a primary system, which is subject to an external excitation, into a secondary’ system, 
which is not a subject to an external excitation. Once the secondary system is excited, the energy flow occurs in 
both directions. 18 The energy flowing back from the secondary to the primary system in conjunction with a viscous 
damping can cause the secondary system response to exhibit an intermittent behavior. 

An aluminum arch with shallow curvature under transverse loading is considered. Comparisons between 
nonlinear modal simulations and numerical simulations in physical DoFs are made to determine the proper modal 
basis selection. Of particular interest is the condition of autoparametric resonance and its effect on the required 
basis selection. 


II. Nonlinear Numerical Simulation in Physical DoFs 

The shallow arch with cylindrical curvature was first analyzed with numerical simulation in physical DoFs to 
give an insight into the characteristic behavior of non-planar structures. In this respect, the simulation is treated as if 
it were an experimental observation. It will therefore help guide the selection of basis functions used in the reduced 
order nonlinear modal simulation, and will ultimately be used to validate those results. A clamped-clamped 
aluminum arch measuring 18-in. x 1-in. x 0.09-in (projected length x width x thickness) and 81.25-in. radius of 
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curvature was considered, as shown in Figure 1. This geometry gave a radius/span ratio of approximately 4.5. The 
following material properties were used: 

E = 10.6 x 10 6 psi, G = 4.0xl0V', P = 2-588 x KT 4 . 

Mass proportional damping was specified corresponding to critical damping of 2.0% for the first symmetric 
transverse mode (at 258 Hz). A random load generation process previously presented 9 was applied, giving a flat 
excitation spectrum from 0 to 1500 Hz. The arch was uniformly loaded along its span in the vertical direction, 
irrespective of the deformation, i.e. follower forces were not utilized. 

The finite element program ABAQUS (version 6.5-1) was used to generate nonlinear displacement time 
histories. The double precision explicit integration scheme with an adaptive time integration step (referred to as 
“element by element” in ABAQUS) was utilized for all analyses. A total time of 2.1384.? was simulated for each 
loading level considered, and the first 0.5? was removed to eliminate the transient response. The ABAQUS model 
consisted of 144 B21 beam elements. All displacement results are presented in the global coordinate system as 
horizontal (u) and vertical (v) components. 

The responses to four different excitation levels yielding slightly to highly nonlinear response characteristics 
were computed. The v and u and displacement power spectral densities (PSDs) were examined at two locations, 
namely at the projected quarter-span of the arch (4.5 in. from the clamped boundary), presented in Figure 2 and 
Figure 3, and at the mid-span (9.0 in. from the clamped boundary), presented in Figure 4 and Figure 5 Displacement 
responses at the quarter-span location (Figure 2 and Figure 3) clearly indicate a significant change in the PSD 
characteristics when passing from the 0.4608 lb/in. to the 0.6517 lb/in. excitation level. In particular, the emergence 
of a peak at approximately 158 Hz is noted. At the mid-span location, the v displacement PSD (Figure 4) does not 
indicate this change in character. Examination of mid-span u displacement PSD (Figure 5), however, indicates 
negligible response amplitude at the 0.2304 lb/in. excitation level, followed by a significant jump in amplitude at 
and above the 0.6517 lb/in. excitation level. The 0.4608 lb/in. excitation level appears to produce a response in a 
transition region. 

The jump in response in Figure 5 was investigated in further detail. Figure 6 presents mid-span u displacement 
time histories ranging from the 0.2304 lb/in. to 1.3033 lb/in. excitation levels. At the excitation level of 0.2304 
lb/in., the mid-span u response is negligible. Previous work indicated this to be the case for a planar isotropic and 
uniformly loaded beam, regardless the excitation level applied. 8 For the loading level of 0.6517 lb/in., the mid-span 


4 



u displacement of the arch exhibits intermittent damped transients. At the excitation level of 1.3033 lb/in., the 
motion is of a permanent nature. A few instantaneous snapshots of a full-field displacement along the span of the 
arch for the 0.6517 lb/in. loading were studied to determine the nature of this motion. Snapshots were chosen from 
the evolved u displacement time history shown in Figure 7, such that the center point displacement was positive (left 
of the neutral position, t = 1.701m), negative (right of the neutral position, t = 4.148 ms) and nearly zero (at the 
neutral position, t = 3.225 ms). Snapshots of the v and u displacements along the arch span for these three conditions 
are presented in Figure 8 and Figure 9, respectively. When the mid-span u response is zero (t = 3.225 ms), the v 
displacement is symmetric and u displacement is anti-symmetric. It is seen in both figures that symmetry of the v 
displacement response and anti-symmetry of the u displacement response is only preserved when the mid-span it 
displacement is zero. When the mid-span u displacement response is positive (t = 1.701m), the v displacement 
response is no longer symmetric and its maximum occurs at approximately the quarter-span point from the right end. 
When the mid-span u displacement response is negative (t = 4.148m), the v displacement response is also no longer 
symmetric, with its maximum at approximately the quarter-span point from the left end. These observations lead to 
the conclusion that both u and v displacements are generally a superposition of both symmetric and anti-symmetric 
behaviors. 

Due to the symmetry of the arch geometry and the loading distribution, it can also be concluded that there exists 
no direct form of excitation that would initiate the anti-symmetric v and/or symmetric u displacement responses, i.e., 
non-symmetric v displacement and/or any mid-span u motion. Therefore this behavior must be initiated by energy 
transfer from the primary system, which is consistent with the geometry and loading. Recalling Section I, such an 
observation fulfills a definition of an autoparametric resonance. The analysis in physical DoFs, however, is unable 
to provide detailed insight as to what modes contribute to the anti-symmetric behavior. Certainly, inspection of 
Figure 2 and Figure 3 can provide some clues regarding the new frequency components introduced with the anti- 
symmetric v displacement response (i.e. at approximately 160, 520, 1070 Flz). Flowever, since the response is 
nonlinear, characterization of this nonlinear behavior in terms of linear modes is appropriate only in the context of a 
nonlinear modal simulation. 

Founded on experience built upon flat and isotropic thin-walled structures analysis, there exists a perception that 
a nonlinear analysis is generally an attempt to mitigate conservativeness or over-design that would result from a 
linear analysis. Flowever, nonlinear curved and/or anisotropic structural analysis can result in a solution that will 
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not be upper-bounded by a linear analysis. 13 ' 14 Figure 10 and Figure 11 present linear and nonlinear vertical root 
mean square (RMS) displacement results at the mid-span and quarter-span locations, respectively. It is clearly seen 
that for the shallow arch analyzed herein, the displacements at both selected locations obtained from the nonlinear 
analysis will exceed the results obtained form the linear calculations. Hence, this structure exhibits a spring 
softening nonlinearity and is non-conservative relative to the linear system. 

III. Reduced Order Numerical Simulation 

The nonlinear reduced order numerical simulation analysis consists of several parts. 8 The linear eigenvectors, 
commonly employed as basis functions, are first obtained from a commercial finite-element program (in this study 
ABAQUS v.6.5-1). Following a transformation of the nonlinear system of equations to modal coordinates, the 
modal stiffness coefficients are evaluated and the resulting coupled system of equations is numerically integrated to 
obtain the modal displacement time history. The modal solution is transformed back to physical coordinates and can 
be post-processed as desired. It should be pointed out that for the purpose of analyzing non-planar structures, the 
formulation previously presented 7,8 ' 11 has not been redeveloped. However, the guidelines for the selection of basis 
functions must be substantially revised. 

A. Modal Coordinate Transformation 

The equations of motion of the nonlinear system in physical DoFs may be written as 

MX(t) + CX(t) + F NL (X(t)) = F(t) (1) 

where M and C are the mass and proportional damping matrices, respectively, and X and F are the displacement 
response and force excitation vectors, respectively. The nonlinear restoring force F NL contains the linear, quadratic 
and cubic stiffness terms. 

A set of coupled modal equations with reduced degrees-of-freedom is first obtained by applying the modal 
coordinate transformation X = Or/ to Equation (1), where q is the modal displacement response vector. The modal 
basis function matrix d> is typically formed from the eigenvectors obtained from Equation (1) using only the linear 
stiffness. Generally, a small set of L basis functions are included giving 

Mq(t) + Cq(t) + F nl ( q t ( t ), q 2 ( t ), ...,q L (t)) = F(t) (2) 

where the tilde superscript represents modal quantities, and 
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F = ® t F . 


( 3 ) 


M = 3> r MO =\l J C = 0> r C<t> = T 2 C r o r J F nl = O t F nl 

B. Indirect Stiffness Evaluation Method 

The previously developed indirect stiffness evaluation method is applied. 7 ' 8 ' 11 To summarize, the nonlinear force 
vector in Equation (2) may be written in the form 

L L L L L L 

T, d 'j1j r = \,2,...,L (4) 

7=1 7=1 k— j 7=1 k-j l-k 

where d , a, and b are the linear, quadratic nonlinear, and cubic nonlinear modal stiffness coefficients, respectively. 
This form reduces the problem of determining the nonlinear stiffness from one in which a large set of simultaneous 
nonlinear equations must be solved to one involving simple algebraic relations. The algebraic relations are obtained 
by solving a series of nonlinear static problems with prescribed displacement fields. 7,8,11 Note that the indirect 
procedure, contrary to the direct methods, 1 " 314 does not neglect in-plane inertia of the system. 

C. Eigenanalysis 

As mentioned briefly in Section III.A, a reduced order method usually employs a subset of selected eigenvectors 
as a set of basis functions to be used in a reduced order stiffness evaluation procedure. As expected and explained in 
Section I, the identification and classification of these eigenvectors for a shallow structure is more involved than for 
planar configurations. It is helpful to categorize the linear eigenvectors that constitute the modal basis by their 
symmetry. For the shallow shell under consideration two categories are defined. Those having symmetric 
transverse displacements, with symmetric rotations and in-plane displacements, are subsequently referred to as ST 
modes. Those having anti-symmetric transverse displacement with symmetric rotations and in-plane displacements, 
are referred to as AT modes. Each eigenvector consists of both transverse and in-plane components, therefore an 
additional criteria is needed in an attempt to recognize whether a particular mode will be contributing more to a 
transverse or in-plane response component. In this work the following classification was adopted: For each mass 
normalized eigenvector, the maximum transverse and in-plane components were determined. If the maximum 
transverse component was greater than the maximum in-plane component, the mode was classified as type T - 
transverse dominated. If the maximum in-plane component was greater than the maximum transverse component, 
the mode was classified as type I - in-plane dominated. Hence, four classes of modes are considered for inclusion in 
the modal basis, i.e., ST-T, ST-I, AT-T, and AT-I. Note that this scheme for classifying mode type is applicable 
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only to the shallow curvatures under consideration here. For deeper curvatures, an alternative classification scheme 
would likely be necessary. It should also be noted that either the full span model, or a semi-span model analyzed 
with both symmetric and anti-symmetric mid-span boundary conditions, are needed to obtain the above suite of 
eigenvectors. 

The identical ABAQUS model was used for both the physical DoFs solution and eigenanalysis for the reduced 
order numerical simulation. Selected modes with their ordering numbers, corresponding natural frequencies, and the 
above-introduced classification are shown in Table 1. For the structure considered, type T modes typically have 
lower natural frequencies than type I modes. Vertical displacement, rotation, and horizontal displacement 
components of the lowest natural frequency eigenvectors of each of the four categories defined above are shown in 
Figure 12 to Figure 15. 

At this point it is also worthwhile to note that the concept of so-called companion 10 or dual 1 112 modes developed 
for a planar structures is no longer needed nor relevant. These bases were previously introduced to obtain 
membrane behavior corresponding to bending modes of the flat structure. Since the modes are linearly coupled for 
curved structures, these degrees of freedom are automatically reflected in the basis. 


D. Modal Stiffness Relationships 

It was previously shown 8 that the linear modal stiffness obtained from the indirect stiffness evaluation method is 
equivalent to that obtained via a standard modal analysis, that is, 

d = K = <& T K<& = \ co; \ (5) 


As expected, there is zero modal coupling through the linear modal stiffness. Additionally, certain relationships 
exist between significant quadratic and cubic nonlinear stiffness terms. 7 ' 8 For the quadratic terms, these are 


a — ^a 2 
2a — a 2 


(6) 


For the cubic terms, these are 


1 1 
^122 

= Kn 

b' 

= b 2 

u 222 

u \22 

K 

= b\ n 


( 7 ) 


Without citing particular values, these relationships are maintained for the present problem. 



An inspection of the nonlinear modal stiffness coefficients can be helpful to reveal couplings that may exist 
between specific types of bases. In the following, it should be recognized that signs of the modal nonlinear stiffness 
coefficients are only affected by signs of the utilized basis functions (in this case eigenvectors), which are not 
unique, and by signs of scaling factors applied to them, which are also arbitrary. Therefore, the absolute value of 
coefficients will be used to identify coupling of modes. Since the quadratic stiffness matrix is three-dimensional and 
the cubic stiffness matrix is four-dimensional, a spot inspection of numerical values may not be particularly 
revealing. Thus, a graphical representation is offer herein. 

Figure 16 presents the absolute values of quadratic stiffness coefficients ci jk for the basis formed out of 4 AT-T, 

4 ST-T, 4 AT-I, and 4 ST-I bases, that is, AT-Ti_ 4i ST-T ,_ 4 AT-I,_ 4 ST-I^. These bases are numbered from 1-16, 
e.g. basis 15 refers to ST-I 3 . The basis selection is discussed in the next section. Each of the 16 subplots has index i 
fixed, and indices j and k are shown on the horizontal and vertical axes of each subplot, respectively. The color bar 
indicates the coupling strength, with dark blue effectively depicting a lack of coupling between the basis functions 
involved. It is seen that similar coupling patterns between three modes are indicated for each row of subplots. For 
example, the top row shows strong coupling between AT-T modes (l</<4) and AT-T modes ( I </<4 ) and ST-I 
modes (13<&<16), and also between AT-T modes (l</<4) and ST-T modes (5</'<8) and AT-I modes (9<k<\2). The 
remaining rows indicate different couplings between three modes. 

As mentioned before, the cubic matrix is four-dimensional, so a two-dimensional graphical representation of 
coefficients b' jkl requires fixing of two indices, say i and j, and varying the remaining two, k and /, on the horizontal 

and vertical axes, respectively. For 16 bases, this procedure would produce 256 two-dimensional plots. Flowever, 
the same trend as previously found for the quadratic matrix was also observed for the cubic matrix, namely, when 
indices i and j each vary within the same type of basis function, the couplings resemble very similar patterns. For 
example, for an i-index variation within the AT-T basis (l</<4), and a y-index variation within the ST-T basis 

(5</<8). plots corresponding to any of coefficients b\ g k ; , are similar. Thus, the cubic modal stiffness coefficients 

can be represented in 16 subplots, as shown in Figure 17. Each subplot indicates different couplings between four 
modes and is representative of sixteen such subplots. 

The specific conclusion to be drawn based on Figure 16 and Figure 17 is that not all possible means of coupling 
are exercised in the particular problem under investigation. At the same time it is observed that none of the basis 


9 



subsets utilized in the analysis is completely uncoupled from the rest of the basis functions. For example, for the 
cubic modal stiffness b \ 2 13 _ 16 13 _ 16 , the coupling between AT-T, AT-I, ST-I, and ST-I modes is effectively zero (see 

b\ lk , in Figure 17). Flowever, a different combination of these same three mode types is coupled by means of the 

cubic modal stiffness b\l 9 _ 12 ,_ 4 (modes ST-I, ST-I, AT-I, and AT-T, see b\l kl in Figure 17). 

In general, it is not possible to determine whether the selected set of basis functions is insufficient from the 
modal coupling plots. However, it is possible to determine which modes would not contribute to the response. If a 
specific mode is completely uncoupled from the rest of the basis, and also is not directly excited, then its presence in 
the basis would not be necessary. The detailed assessment of the basis set sufficiency is reserved for the following 
section. This assessment is based on reconstructed physical displacements following computation of the modal 
response, in a process which involves modal summation. 

E. Reduced Order Analysis Results 

For the shallow curvature considered, two factors affected the initial modal basis selection. The first was earlier 
experience with symmetric planar structures under symmetric transverse loadings where the initial basis consisted of 
low frequency bending modes. 7 " 12 For the curved configuration, an analogous basis consisting of transverse 
component dominated (type T) modes was chosen. Like the earlier experience, the addition of higher frequency in- 
plane dominated modes (type I) was expected to improve the prediction of peak magnitudes and their broadening in 
a highly nonlinear regime. 8 The second selection factor was based on observations made through the simulation in 
physical DoFs. There, the need for both ST-T and AT-T modes was indicated for autoparametric resonance 
conditions. Thus the modal bases chosen consisted of all four types of modes. In the following, an abbreviated 
convention describing the set of modes applied in a particular analysis was used. The convention consists of a four- 
digit sequence containing the number of included ST-T, ST-I, AT-T, and AT-I modes, in that order. Note that this 
ordering is not based on modal frequency, and hence differs from that of Figure 16 and Figure 17. For example, 
4044 indicates a basis consisting of four ST-T, zero ST-I, four AT-T, and four AT-I modes. 

Two excitation levels were selected for investigation; 0.2304 lb/in. was shown to give rise to a slightly nonlinear 
response, and 0.6517 lb/in. produced a highly nonlinear response. As illustrated in Figure 6, the lower excitation 
level does not produce an autoparametric response, while the higher level produces an intermittent autoparametric 
behavior. For both loading levels, the following suite of basis functions was considered. The 4040 basis was 
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selected as the baseline as it contained only type T modes, i.e. four ST-T and four AT-T modes. The 4440 basis was 
used to determine the effect of ST-I modes on the baseline. Similarly, the 4044 basis was used to determine the 
effect of AT-I modes on the baseline. The effect of both types of in-plane dominated modes on the baseline was 
studied using the 4444 basis. Two additional bases were considered to investigate the contribution of only ST-T 
modes in the 4000 basis, and the contribution of only AT-T modes, in the 0040 basis. These last two bases were 
ideal for investigating the response to direct excitation. The ability to interrogate the nonlinear response in this 
manner is a unique capability of the nonlinear modal simulation. An analysis in physical DoFs can only 
approximate this ability by limiting the frequency bandwidth of the excitation field. The effect of nonlinear modal 
coupling however remains in that case. 

Having the nonlinear force vector in Equation (4) fully defined with any set of the basis functions described 
above, the coupled modal nonlinear equations of motion given by Equation (2) were numerically integrated using a 
fourth-order Runge-Kutta method. 19 Consistent with the analysis in physical DoFs, a total time of 2.1384.? was 
simulated and the initial 0.5? of transient response was removed. A fixed time integration step of 2 /Js was used. 
The resulting modal displacement time histories were transformed back to physical coordinates using the inverse 
modal transformation. 

Horizontal and vertical displacement PSDs obtained via reduced order analysis were compared with those 
obtained via the analysis in physical DoFs at two locations along the span of the arch, namely at the center and 
quarter-span nodes. The need to select more than one validation point was based on the observation from the 
analysis in physical DoFs that the change in nonlinear response manifests itself in different ways at different 
locations. At the quarter span location, the emergence of an additional peak was noted in Figure 2 and Figure 3. At 
the mid-span node, the u displacement component exhibited a sudden jump in response level with the onset of 
autoparametric resonance. 

Figure 18 to Figure 20 present the comparison of the reduced order solutions obtained with different bases and 
the solution in physical DoFs at the 0.2304 lb/in. excitation level. For these and subsequent plots, significant peaks 
are annotated with their corresponding modal component to aid the discussion. It should be understood that these 
annotations do not imply the presence of linear modes in the nonlinear response. The u displacement at the mid- 
span is not shown, since all analyses yielded numerical zero. All reduced order displacement PSDs compare well 
with the physical DoFs solution and demonstrate behaviors due solely to ST modes. This is evidenced by the fact 
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that the inclusion of AT-T and AT-I modes in the basis do not alter the computed response, e.g. compare 4000 with 
4040 and 4044. In fact, the 0040 basis resulted in a numerically zero response at both locations for both response 
components, indicating that the AT modes are not directly excited by the externally applied loading. The only 
noticeable difference between modal bases is apparent at the quarter-span location in the first sub-harmonic of the 
lowest ST-T mode (133 Hz). Upon close inspection, it is seen that the behavior is represented slightly better with 
the 44xx bases (4444 and 4440), than it is with any of the 40xx bases (4044, 4040, and 4000). Therefore, even for 
this slightly nonlinear response regime, the inclusion of ST-I modes is shown to improve the behavior related to the 
ST-T modal response. As evidenced by Figure 16 and Figure 17 introduced in Section III.D, this is due to the 
presence of cross coupling of nonlinear stiffness terms, in this case between the ST-T and ST-I mode pairs (see 
a 5 , k K in Figure 16, and Z^^and b\ bkl in Figure 17, for example). Note that if the excitation bandwidth were 

sufficiently increased, the ST-I modes would also be directly excited. From this load level, we conclude that the ST- 
T modes (though direct excitation) and ST-I modes (through coupling with ST-T modes) constitute the primary 
system. 

Figure 21 to Figure 28 show the results obtained at an excitation level of 0.6517 lb/in. For clarity, each set of 
solutions is displayed in two separate figures. One figure indicates the physical DoFs, 4444, and 4044 solutions. 
The second figure shows the physical DoFs, 4440, 4040, and 4000 solutions. 

At the mid-span location, the v displacement response (Figure 21 and Figure 22) is dominated by the peak 
originating from the first ST-T mode. The response is best modeled by the 4444 and 4440 bases. Since the PSD 
response does not indicate a presence of peaks originating from anything other than ST-T modes, this result is found 
consistent with the conclusions based on 0.2304 lb/in. case, where the comparison of 44xx bases with the physical 
DoFs solution was found superior to the 40xx bases. The 40xx bases (4044, 4040, and 4000) over-predict the 
degree of nonlinearity. This is evidenced through excessive stiffening (peaks moved to the higher frequencies 
compared to physical DoFs) of the first, third and fourth ST-T-related peaks, and excessive broadening, especially of 
the fourth peak. The second ST-T-related behavior is not significantly present at this location. The solution for the 
0040 basis was numerically zero, indicating that the AT-T modes were not directly excited. 

Mid-span u displacement results (Figure 23 and Figure 24) offer a different perspective. The response exhibits 
peaks originating from the AT-T modes. The 0040 basis yielded a numerically zero response, again indicating that 
the AT-T modes are not directly excited by the externally applied loading. The excitation provided to the AT-T 
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modes is therefore applied indirectly through autoparametric resonance. Thus, the AT-T modes (through 
autoparametric resonance) and the AT -I modes (through coupling with the AT-T modes) constitute the secondary 
system. The modal bases containing ST and both AT-T and AT-I modes are therefore expected to correspond best 
with the physical DoFs solution for this condition. The solution for basis 4044, shown in Figure 23, corresponds 
fairly well with the physical DoFs solution. Flowever, the magnitude of the first AT-T peak is too low, and it 
demonstrates over-stiffening. Further, the comparison beyond approximately 750 FIz frequency range is deficient. 
The addition of ST-I modes in basis 4444 activates the ST-T to ST-I coupling and improves the comparison 
throughout the frequency range. Bases which do not include AT-I modes (4440 and 4040) suffer as the AT-T to 
AT-I coupling is not represented, see Figure 24. The 4440 basis captures the behavior of the first AT-T peak better 
than the 4040 basis, likely due to the ST-T to ST-I coupling. Beyond about 300 Hz, both 4440 and 4040 bases 
compare poorly, as the response level generally under-predicts the physical DoFs solution except in the frequency 
range between 1050 and 1200 Hz, where the contribution of the third AT-T is exaggerated. It can be concluded that 
the higher frequency modes of both kinds, ST-I and AT-I, are needed for an accurate mid-span horizontal 
displacement prediction when the structure undergoes autoparametric resonance. However, inclusion of the AT-I 
modes is more important than inclusion of the ST-I modes, as evidenced by the improved comparison of the 4044 
basis relative to the comparison of the 4440 basis. This observation further substantiates the importance of the 
coupling between the ST-T and ST-I modes (as indicated above), and between the AT-T and AT-I modes (see 

in Figure 16, and b l l2kI and tf 6tl in Figure 17, for example). Note that basis 4000 resulted in an essentially zero 

response, indicating that the basis is incapable of representing the secondary system and capturing the 
autoparametric interactions. 

At the quarter-span location, the v displacement component (Figure 25 and Figure 26) is best captured by 4444 
basis. Basis 4000 does not have means of representing the anti-symmetric behavior, and consequently provides a 
solution that fails to capture the behavior of the secondary system. This is particularly apparent in the lack of the 
peak originating from the first AT-T mode, including its sub-harmonic peak. Among those bases that are capable of 
reflecting the anti-symmetric behavior, basis 4040 is of the lowest quality. It poorly represents the peak originating 
from the first AT-T mode in that it is over-stiff (its peak magnitude is too small and is shifted toward higher 
frequency). Additionally, its AT-T sub-harmonic is too small while the contribution of the third AT-T mode is 
exaggerated. The quality of solutions from bases 4440 and 4044 consistently falls between the best (4444) and the 
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worst (4040). Overall, basis 4440 matches the peak magnitudes and broadening better than the 4044 basis. The 
exception is in the frequency range related to the third AT-T mode (approximately 1050 - 1150 Hz), in which the 
solution from the 4440 basis is vastly exaggerated, while the solution from the 4044 basis over-predicts this peak 
only slightly. 

At the quarter-span, the u displacement component (Figure 27 and Figure 28) obtained with 4444 basis is again 
superior to all the others. It very accurately matches peaks originating from both the ST-T and AT-T modes. It is 
particularly worthwhile to note that the 4444 basis is the only one that does not exhibit over-stiffening of the first 
AT-T related peak, and accurately captures the second and the third AT-T related peaks. Remarks made regarding 
the performance of 4440, 4040, and 4000 solutions for the vertical displacement component apply to the horizontal 
displacement component as well. The only exception is that the 4044 basis under-predicts the peak originating from 
the third AT -T mode for the u displacement component. 

IV. Conclusions 

The effect of modal basis selection on the displacement response predicted using a nonlinear modal simulation 
was investigated. For the shallow arch considered, it was found that: 

> At a sufficiently high excitation level autoparametric resonance can exist; 

> For a symmetric curved structure and symmetric loading distribution, the primary system consists of ST-T 
and ST-I modes and these modes may be directly excited by external loadings with sufficient bandwidth; 

> For a symmetric curved structure and symmetric loading, modes AT-T and AT-I are never directly excited 
(regardless of the excitation level and bandwidth); these modes constitute the secondary system and are only 
indirectly excited through the nonlinear modal interactions at a sufficiently high response regime; 

> For non-planar symmetric structures which may exhibit autoparametric resonance, the inclusion of both ST- 
T and AT-T modes is the necessary minimum to represent directly and indirectly loaded components; 

> Higher frequency mode variants (ST-I and AT-I) provide a significant enhancement of the solution quality 
in the highly nonlinear/autoparametric regime; for all the cases studied at the high excitation level, inclusion 
of all mode variants (ST-T, ST-I, AT-T, and AT-I) always provided the solution that compared best with the 
physical DoFs analysis; 
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> The inclusion of all mode variants at any location or excitation level does not adversely affect the response 
prediction; when mode identification is difficult, it should be always “safe” to expand the set of basis 
functions; 

> In the highly nonlinear/autoparametric regime, modes ST-I appear to enhance ST-T-originating response 
features, while AT-I appear to enhance AT-T-originating response features; therefore, modal pairs: ST-T 
with ST-I, and AT-T with AT-I exhibit a strong coupling; 

> Validation of a reduced order solution versus physical DoFs analysis and/or experimental results should be 
performed at more than one location on the structure and/or by inspection of more than one response 
component; e.g., if only transverse deflection at the center point would be considered, basis 4000 might have 
appeared sufficient; 

> Reliable evaluation of the quality of the solution as a function of modal truncation is difficult during the 
intermediate steps of the analysis in modal coordinates (i.e. based on modal stiffness coefficients); an 
accurate assessment can only be provided following the modal response computation and inverse 
transformation to physical DoFs; 

> Nonlinear analysis of shallow curvature structures may yield non-conservative results when compared with 
linear analysis. 

> The ability to interrogate the nonlinear response by altering the modal basis is a unique capability of the 
nonlinear modal simulation; analyses in physical DoFs will always reflect nonlinear modal coupling. 

The results presented in this paper were for a particular shallow curvature. A possible area for future work is to 
determine the onset of autoparametric resonance as a function of curvature. A related interest is to extend the 
assessment of modal basis selection to structures having greater curvature and to those having linear coupling due to 
material anisotropy. At a minimum, it is expected that the present modal classification scheme would need 
modification for these cases. 
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Figure 1: Geometry of shallow arch. 
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v Displacement PSD, in 2 /Hz 



Figure 2: Vertical displacement PSD at quarter-span location. 
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u Displacement PSD, in 2 /Hz 



Figure 3: Horizontal displacement PSD at quarter-span location. 
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Figure 4: Vertical displacement PSD at mid-span location. 
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Figure 5: Horizontal displacement PSD at mid-span location. 
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Figure 6: Mid-span horizontal displacement response time histories. 
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Figure 7: Mid-span horizontal displacement time history at 0.6517 lb/in. distributed loading. 
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v Displacement, in. 



Figure 8: Vertical displacement field for selected time instances at 0.6517 Ib/in. distributed loading. 
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Figure 9: Horizontal displacement field for selected time instances at 0.6517 Ib/in. distributed loading. 
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Figure 12: The first type T anti-symmetric transverse mode AT-Ti (158Hz). 
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Figure 13: The first type T symmetric transverse mode ST-Ti (258 Hz). 
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Figure 14: The first type I anti-symmetric transverse mode AT-Ij (5.6 kHz). 
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Figure 15: The first type I symmetric transverse mode ST-Ix (11.2 kHz). 
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Figure 16: Quadratic nonlinear modal stiffness coefficients. 
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Figure 17: Cubic nonlinear modal stiffness coefficients. 
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Figure 18: Mid-span vertical displacement at 0.2304 lb/in. distributed loading. 
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Figure 19: Quarter-span vertical displacement at 0.2304 ib/in. distributed loading. 
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Figure 20: Quarter-span horizontal displacement at 0.2304 lb/in. distributed loading. 
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Figure 21: Mid-span vertical displacement at 0.6517 lb/in. distributed loading. 
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Figure 22: Mid-span vertical displacement at 0.6517 ib/in. distributed loading. 
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Figure 23: Mid-span horizontal displacement at 0.6517 lb/in. distributed loading. 
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Figure 24: Mid-span horizontal displacement at 0.6517 lb/in. distributed loading. 
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Figure 25: Quarter-span vertical displacement at 0.6517 ib/in. distributed loading. 
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Figure 26: Quarter-span vertical displacement at 0.6517 Ib/in. distributed loading. 
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Figure 27: Quarter-span horizontal displacement at 0.6517 lb/in. distributed loading. 
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Figure 28: Quarter-span horizontal displacement at 0.6517 lb/in. distributed loading. 
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Table 1: Natural modes - classification and frequencies. 


Mode 

Designator 

Mode 

Number 

Frequency 

(Hz) 

ST-T, 

2 

258.17 

st-t 2 

3 

400.41 

st-t 3 

5 

773.72 

st-t 4 

7 

1,426.0 

ST-T, 

22 

11,225 

st-i 2 

33 

22,437 

st-i 3 

42 

33,638 

ST-I 4 

51 

44,826 

AT-Tj 

1 

158.25 

at-t 2 

4 

513.19 

AT-T 3 

6 

1070.1 

AT-T 4 

8 

1828.0 

AT-Ij 

15 

5,624.3 

at-i 2 

28 

16,831 

AT-I 3 

38 

28,039 

at-i 4 

47 

39,235 
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